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Form  of  the  Problem-  Using  matrix  notation,  the  standard  problem 
may  be  written 

(1)  Minimize  wb  *■ 

subject  to  wA°  >  c° 

and  w  *  0=  w  integer 

where  b  is  an  m  X  1  eolmn  vector,  is  a  scalar,  e°  and.  0  are  IXn  row 
vectors,  A°  is  an  mXn  rectangular  matrix,  and  w  is  a  lXn  row  vector 
whose  values  wn  wish  to  find  so  as  tc  achieve  the  minimisation  objective 
and  satisfy  the  constraining  conditions  of  (1),  Th«  fact  that  w  must  be 
lnt©gB*l  distinguishes  the  problem  from  the  general  linear  programing 
problem  which  allows  w  to  have  fractional  components . 

For  the  purpose  of  this  paper  we  define  the  mX(m  ■»  n)  augment ed 
matrix.  A  and  the  II (m  *  n)  augmented  vector  «  by 
A  »  (A°  I),  e  *  h°  0), 

(where  .1  Is  the  mXm  identity  matrix)  and  rewrite  (1)  as 

(2)  Minimize  wb  t 

subject  to  wA  >  ev  w  integer. 

We  assume  without  loss  of  generality  that  the  augmented  matrix  (b  A)  is 
lexi cograph icaliy  negative  by  row*  for  if  K  is  not  it  may  readily  be 
made  6o  {sm  (ij)..  Following  Gomory's  terminology  [3],  we  will  call 
our  method  an  all-integer  algorithm,  for  we  additionally  require  that  all 
elements  of  A,  b„  and  <c  be  integral,  or  at  least  conmensurable ,  In  practical 
terms  this  is  of  course*  no  restriction  s.in?e  numbers  are  represented  in 
the  computer  with  finite  decimal  expansions  in  any  case. 

Tools  of  the  Algorithm , 

1  A  set  of  transformations  which  will  change  the  problem  into  a 
new  problem  in  non.negatt.ve  variables  so  that  any  optimal  integer  solution 


to  th«  new  problem  provide®  an  optimal  integer  Solution  to  th«  original,, 
and  •^onvenuly,  We  will  wall  these  trans  formats ons  elemental  transformation 

2.  A  prooedaral  rule  (eo spied  with  a  rule  of  choice)  for  applying 
the  elemental  transformations  in  order  to  create  a  new  problem  containing 
a  suhaatrix  of  a  special  form,  which  we  will  call  the  bounding  form,, 

3«  A  process  called  the  bound  escalation  method  for  operating 
on  the  bounding  form  to  supply  lower  bound  values  for  some  set  of  the 
problem  wariafeSes., 

Thus  the  algorithm  may  bo  roughly  sketched  as  follows. 

1.  Apply  a  series  of  elemental  transformations  to  obtain  an 
equivalent  problem  in  new  variables  which  exhibits  a  bounding  form  .matrix  B. 

2.  Apply  tha  bound  escalation  method  to  B„  At  the  end  of  the 
process  the  lower  bound  values  assigned  to  a  subset,  of  the  problem  variables 
will  satisfy  all  the  constraints  associated  with  B. 

'$*  Adjust  the  s  vector  to  reflect  the  assignment  of  lower  bounds 
established  in  2,  If  a  becomes  nonpositive,,  the  problem  is  solved, 
Otherwises,,  return  to  1  and  repeat. 

Several  features  of  the  algorithm  may  be  noted.  First*  th#  problem 
is  sol wed  directly,  i.e.,  no  reference  is  made  to  the  dual.  Second,  there 
is  no  pivoting  proceed  In  its  customary  form. 

Thus,  the  elemental  transformations  are  applied  until  the  problem  .Vs 
ready  for  the  bound  escalation  method,  and  then  the  machinery  for  the 
latter  is  set  into  action-,  Those  who  wish  may  relate  these  two  steps  to 
a  fo«u  of  deferred  pivoting  and  abbreviated  pivoting*  respectively,  bv.t 
attempts  to  salvage  the  pivoting  concept  are  inessential  to  understanding 
the  process.  Third.  .no  ua«  is  made  of  slack  variables  to  transform 
ine«f\«siione  into  equations.  Fourth,  the  netiio-)  v.  -rb,  guns rally  j.o 


satisfy  some  set  o'*  constraints  simultaneously  with  the  bound  escalation 
method  .rather  than  talking  the  mow  narrow  lamedtate  view  of  satisfying 
a  single  constraint,.  .Fifth.  the  bound  escalation  method  leases  all 
constraints  of  the  bounding  form  satisfied,  whereas  the  pivoting  process 
of  other  integer  algorithms  may  not-  in  one  stop  completely  satisfy 
the  constraint  to  which  they  are.  applied.  Sixth,  because  there  is  no 

e 

customary  pivot  operation,  a  choice  among  eligible  pivotal  constraint 
is  replaced  by  a  choice  of  another  sort,  i.e»6  that,  of  the  sequence 
of  elemental  transformations  with  which  to  establish  a  bounding  form. 

We  will  now  lay  down  the  basis  of  the  algorithm.  Proofs  of  the 
lemmas  to  follow  will  be  found  in  appendix  I„ 


l .  ?hs  Elemental  Transformation 3 . 

Consider  the  set  of  transforms  T  «  {  T^8,  T^8  },  rys  *  1,  2,  > ns, 


„ 

*  s.,  where  w«  define  the  components  t# ,  of  the  mXm  matrix  T.  by 

1.J  K 

(  -1  if  i  ■»  r}  j  *  s 


■‘S'  f 

•>  J 


|  6^  ^  (the  Krone<tker  delta)  otherwise. 


For  example,,  if  b  »  4  w  may  write 
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t2* 


it  !t  littecUately  seen  that  the  matrix 


0 
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°\ 

-1  | 

c  J 

! 
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m  *!.*  «s  .t  except  tr,  rwv.  r  ?i>  wfcl-h  uass  the  components  of  i; 

. .  .  0 

ar«  by 

.Jr 

,•  -  ■ m  -v  n 


r -  ■  a  j  * 

«  l«a*:n‘  s  ..•»*  T  a.~»  s.i  «..»»<  a)  t  rang  forma*  lone,  and  rcneist  simply 

of  rrw  e,Jd‘ t'o.na  and  mi'  ration  .  in  the  A  matrix ,  as  ran  be  seen  by  tn« 


4 


foregoing  remarks , 

We  now  introduce  the  following  two  problems 

(3)  Minimize  wb  * 

subject  to  wA  >  e 
and 

(4)  Minimize  z(Rb)  * 

subject  to  z(RA)  >  », 

which  we  relate  by  the  following  two  lemmas „ 

Lemma  .1,  Assume  that  (3)  and  (4)  have  finite  optima,  and  that  R  may  be 
decomposed  into  a  product  of  elemental  transformations,,  Then  the  vector 
11*  is  a  feasible  (optimally^ solution  of  (3)  if  and  only  if  a*  *•  w*R  'x 
is  a  feasible  (optimal)  integer  solution  of  (4).-, 

Lemma  2,1  Let  R  be  any  transformation  for  which  Lemma  1  is  universally 
true  for  the  standard  integer  programming  problem  of  formulation  (1) 

for  any  A°P  b,  and  s°  satis.fjd.ng  the  finite  optima  .restriction), 
and  let  Rq  be  obtained  from  R  by  re indexing  two  columns  of  R,  Then  either 
R  or  Ry  can  be  expressed  as  a  product  of  elemental  transformations. 

Since  reindexing  two  columns  of  R  corresponds  to  re.index.ing  two 
rows  of  the  augmented  matrix  (b  A)„  and  since  the  latter  merely 
•changes  the  order  in  which  we  list  the  variables  without  changing  the  basic 
problem  in  (3)»  Lemma  2  shows  that  there  is  a  form  of  universality 
inherent  in  the  elemental,  transformations  in  their  application  to  the 
Integer  programming  problem,  • 

For  the  present  algorithm  we  wish  to  restrict  R  somewhat  mor«  than 
in  Leiisas  1  and  2„  To  facilitate  the  ensuing  discussion  we  record  problems 
(3)  and  ( 4 ) .; .truth*  respectively,  in  the* tabular  forme 


*> 


iv) 


- — 

-b 

A 

bo 

e 

U" 


and  define  the  augmented  matrices  corresponding  to  the  upper  portions 
of  (3")  and  (4')  to  be  the  tableau  matrices  of  problems  (3)  ana  (4). 

We  observa  that  the  first  column  of  the  table  (3a)  corresponds  to  the 
objective  function  wfc  »  b^»  and  that  the  successive  columns  identify 
the  constraints  given  by  the  matrix  inequation  wA  >  «?:„ 

One  of  the  restrictions  we  wish  to  put  on  R  is  that  the  tableau 
matrix  of  (4)  is  1 exi  «o graph ic  al ly  negative  (by  row)  whenever  the 
tableau  matrix  of  (3)  is  lexicographically  negative.  The  other  is  that 
the  necessity  for  any  feasible  solution  of  (3)  to  be  nonnegative  implies 
that  fifty  feasible  solution  of  (4)  will  be  also. 

If  we  factor  R  into  alamenfca!  transformations,  we  see  that 

we  are  assured  of  both  cf  these  provided  each  successive  factor  of  R 
assures  them.  For  after  each  step  we  may  redefine  A  and  b  to  equal 
the  matrix  and  vector  RA  and  Rb  just  obtained,  and  re apply  the  result  , 

By  means  of  this  masoning  we  take  feare  of  the  nonnagativity  restriction  In 
Lemma  3„  Assume  that,  any  feasible  solution  of  (3)  must  be  nonnegatlve, 
a&d  let  z*  denbbe  any  feasible  solution  of  (4).>  Then  z*  >  0  is  implied 
by  either  of  the  following, 


b 


t i )  K  *  rjS  for  any  jr  and  e  fr  *  s  )•• 

(it)  R  «  and  fr.hena  exists  a  j  such  that 

c.  >  0„  a.  ,  <  0  for  1  4  Sr,  and  --a  .  >  a  .  *>  0, 

J  ”  ij  ~  rj  -  sj 

Lemma  3  says  that,  we  may  always  assure  that  z*  >  0  by 
subt raet log  one  row  from  another  jjn  the  tableau  matrix.  If  we  add  one 
row  to  another  we  may  still  be  assured  that  z*  >  0  provided  we  identify 

a  j  suffih  that  (a)  Cj  >  0,  and  (b)  the  jth  column  of  A  will  have  all 

components  nonpositive  after  the  addition  except  ag^. 

This  special  form  of  the  jth  column  of  A  in  efon junction  with  >  0 
is  of  particular  interest  to  us.  It  provides  the  fundamental  unit  of 
the  framework  on  which  we  operate  with  the  algorithm  to  converge  to  an 
optimal  solution.  Drawing  on  Lemma  3  we  now  show  how  we  may  manufacture 
this  type  of  column  by  a  sequence  ajf  elemental  transforms,  simultaneously 
preserving  the  desired  restrictions  on  the  problem  form. 

Lemma  Assume  that  (3'  has  finite  optima,  satisfies  the  monnegatlvity 
and  lexicographic  ordering  restrictions,  and  that  the  component  Cj  of  e 
is  positive.  Further  assume  that,  the  tableau  matrix  consists  entirely 
of  integers.  Then  the  following  method  defines  an  R  in  a  finite  number 
of  steps  so  that  problem  U)  satisfies  the  same  restrictions  and  so  that 
the  Jth  column  of  HA  xntalna  exactly  one  positive  component. 

1,  Begin  with  R  ->  f 

2  Select  a  positive  component  of  (a,,),  whew  Iasi.  (a,,)  denotes 
the  Jib  rchsmo  of  A 

If  there  are  n-  other  positive  elements  in  (a  . )  the  procedure 
is  completed,.  Otherwise  pick  a  sedond  positive  component  of  (a,  .}  and 


define  the  subscripts  r  and  »  so  that  a  ,  Is  t.he  component,  associated 

r,  i 

with  the  iexicographt tally  smaller  row  of  the  tableau  matrix  and  a  .  is 


V 

the  component  associated  with  the  larger  row „ 

4»  Redefine  A  and  b  to  be  Tf®  A  and  T^s  b  ( .1  „  e  ,  subtract 
row  s  from  row  r  in  the  tableau  matrix),  redefine  R  to  be  T*8  R„  and 
return  to  instruction  2, 

it  Is  evident  that  therw  may  be  a  variety  of  ways  for  reducing  the 
column  { a , j  '•  to  the  desired  form  by  the  method  of  lemma  4,  some  of 
will  c?h  may  be  more  efficient  in  terms  of  the  number  of  steps  requi  red 
than  others.  One  limned  late  way  that  would  generally  reduce  the  number  of 
steps  would  be  to  replace  in  .Instruction  4  by  (T^  )  Q  where  h  Is  the 
largest  integer  multiple  of  row  »  which  when  subtracted  from  row  t 
will  leave  the  resulting  tow  vector  lexicographically  negatise,  However,-, 
ae  we  shall  see,  the  different  ways  of  reducing  (a,.  , ,  also  may  be  more 
or  } was  efficient  in  terms  of  the  extent  to  which  we  can  exploit.  th» 
structure  of  the  resulting  tableau  matrix.  Hence  at.  this  point  w<* 

■choose  not  to  be  restrictive, 

■  * ••  The  Bounding  form  and  the  Bound  Escalation  Reined, 

Let  1)  bo  a  matrix  whose  columns  correspond  to  some  subset  of  the 
columns  of  Ar  and  let  d  be  the  row  vector  which  corresponds  to  -  in 
exactly  the  same  way  that  Sxes  D  corresponds  to  A  Further  suppose  {i ) 
each  c k;nn  of  D  contains  exactly  one  positive  component and  (H  at  least 
one  of  Uifc  entries  of  d  ’.s  positive.  Finally,  lec  B  be  ‘Vie  matr'y. 
obtained  from  0  by  eliminating  all  rows  in  which  no  positive  .*lment 
appears  Then  we  define  B  to  be  a  bounding  fora  of  A  - 

fr,  t!;p  extreme,  H  may  consist  of  >x  single  -ol-'mu  of  ft  and  ?•  a  single 
l»8h:  vt?  or.nxme.nt  of  the  s  vector,  which  Is  the  "onfigura’.ion  which 

4> 

..-;mr,a  .,  snows  bo*  fo  manufacture  .  In  this  rase  the  ma-i  ri.x  B  -ons  sts 
oVs  ’.he  single  positive  component,  of  0,  We  observe  that  the  inequeUou 


wi;  r  ■.-o;-' ?  >  ;>hiiply  of  «  «  ui  set  of  ♦,|m<  n  .».*&  -.nS  s  defined  ; 

«A  >  u  ,  Ths  bound  ss-ralf^ioii  pro t rdur*  >s  brusssl  .-n  the  rrirripi*  i  n&t 
we  nay  shrink  v  to  x  in  the  same  way  that  D  was  reduced  to  8»  and  irvspia ,■ -r 
wO  >  d  by  :cB  j_  dP  so  that.,  whenever  we  find  an  x*  to  satisfy  the 
latter,  we  have  inpliaitly  a  w*  whifh  satisfies  the  former.  Moreover., 
the  form  of  x  13  enables  us  to  find  an  x*,  hence  a  w* <  suvh  that  the 
w«t. raining  relation  w  >  w*  must  be  satisfied  by  any  feasible 
solution  of  (?}„  'die  following  1  across  show  how  sut?h  a  relation  (which  we 
wi.ll  soon  show  how  to  exploit)  may  be  developed,. 

Lemma  t.  Let.  l)t  8.  d,  and  x  he  given  as  above  We  assume  that  (?)  has 
finite  optima  with  all  variables  constrained  nonnegative  Then  l«  n 
finite  number  of  steps  the  following  procedure  will  fine  fens trained 
lower  bound  values  for  the  roaponenM  of  «  which  vrlJi  satisfy  wC  >  d- 
L,  Let.  x  *  0 

2.  If  all  components  of  d  are  nonposit’ve,  go  to  instruct !on  (, . 

Otherwise,,  pldc  a  positive  component  <>  say  (For  wcplkUnesa,  a 

rcascnaftle  rule  is  to  let  d  ^  be  the  largest  posit ' ve  component, ! 

?..  fmTemant  x^  by  xsM^.  <  d^/b^^-t  where*  b^ ,  >s  the  mSque 

poslfctM  element  it.  the  ,1th  voltaon  of  B,  Kedeflne  1  to  be 

i  *  d  ,/tv  >  (b,  where  (b,  ,)  is  the  kuh  row  of  £  and  return  to  2 
J  £.,1  *'  K 

4..  'ilia  lower  bounds  for  the  components  of  w  wh .<»>  c-.ir^spoml  to 
comrouanfes  of  x  are  given  by  ths  x  vector,  uw  r«aa ’ring  l e uv  nds 
being  0.  iAi  least  one  of  the  components  c-f  x  r. « 


(The  Next.  Page  Js  Page  10- ) 


V'.  e  iviJ  1  illustrate  {.he  metfoci  of  !  emma  ?  w  i  !;  the 


fol  owing  example  problem,  it 1  ready  in 


tabular  for.Dc. 
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front  (  e  lass  '  wo  columns  we  identify  a  bounding  form  w:  ich 
we  $e.  up  in  a  smaller  (able  to  demonstrate  t,he  method, 


•  1 

24j 

3 

-U, 

I  10  -si 

-2 

52 

1 

-20 

-2 

-5 

•he  successive  adjustments  of  1 ->e  <i  vector  are  shown  in  the 
?■  eld  \  *•  ion « 1  rows  below -the  bounding  form  and  the  original  d 
vec: or-  Beside  each  vector  is  che  i nc remen •  of  *  he  varia- 
bid  w!  ich  created  that  vector  out  of  he  Previous  one.. 

Hence  we  end  up  vr;th  vv1  =  3  and  w-.  =  .5.,  Ir  may  be  veri¬ 
fied.  by  substitution  that  these  va  ’u>*s,  sat  isfy  tne  ccnsi  faints 
assuci-ied  wi;'h  f  he  bound  ing’"fonRf  and  in  fact  in  this  case 


«-  S ♦  - 

all  ;  y«-  cons'  ca  in  t  s  of 

the  problem.. 

ft  : 

i  >  o  ‘he  : r t  c on ure  . i u s i 

i 

,y  i  v  .n  is  s u f j c  i  t- 1 1 1 

• 

!o  find 

i 

/ 

\ 

/ 

'he  desired  lower  bound  values  of  w,  the  complete  bound 
escalation  method  is  designed  to  shortcut,  this  procedure 
by  exploiting  certain  properties  of  the  bounding  form. 

We  tum  to  a  considera  ion  of  these  nronerties  with  the 
following  definitions, 

LeS  8  and  E  be  bounding  forms  of  h  S‘>ch  thaf,  B  is  a 
submatrix  of  E  (or  the  same  as  £)„  Then  we  will  call  B  a 
subform  Of  £, 

If  each  row  of  the  bounding  form  E  has  exactly  one 
positive  clement,  we  will  call  E  a  prime  bounding  form,, 
Similarly,  B  will  be  called,  a  prime  subf orm  of  any  bounding 
form  £  if  it  is  a  subform  of  £  and  a  prime  bounding  forme 

We  note  that  ever'’  bounding  form  has  at  least  one  nrime 
subform,  and  also  'ha*  every  subform  of  a  nrime  bounding 
form  is  a  prime  subfonn. 

In  fie  following  examples  of  matrices  with  their  as¬ 
sociated  d  vectors  below  them,  (a)  defines  a  bounding  f orm0 
(b)  defines  a  prime  subform  of  (a)f  and  (c)»  though  it 
may  be  abstracted  from  Ca)s  does  nof  define  a  bounding  form 
a?  all  since  none  of  the  components  of  its  d  vector  are  no s» 
i'ive.  Fin?Il«s  (d)  fails  to  define  a  bounding  form  on  two 
counts;  rhe  firs1  column  contains  more  than  one  positive 
component  and  he  second  row  con  ains  none.  However,,  tie 
form  of  the  las  two  columns  is  such  that  we  may  permissibly 
create  «.  bounding  form  out  of  ’hem  by  removing  the  second  row. 
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Lemma  6.,  Let  ii  be  a  prime  bounding  form,  and  B  a  subform  of 

B„  Let  d  be  associated  ':;if L  B  as  before*  and  define  the 

vector  r  =  oBvl  i hen  there  exists  a  unique  sub  form  B*  of 

*} 

&  sue);  that  (i)  consists  on  J  y  of  nonnegative  components* 

(ii)  each  Xy*  is  positive,  (iii)  if  w*  is  any  feasible 
solution  of  (3),  then  where  the  index 

corresponds  to  fc  as  the  indices  of  c  correspond  to  those  of 
d*,-,  (iv)  if  any  other  subform  B  of  ii  satisfies  properties  (i) 
and  (iii) v  then  B*  implies  a  value  for  each  of  the  components 
of  w+  (in  the  manner  of  (iii))ihat  is  at  least  as  large 
as  implied  by  B 

Lemma  7  „  We  use  •'he  not  a '•ion  of  the  preceding  lemma,  and 

let  h  be  the  vector  associated  with  E  as  d  is  associated 
with  30  Then  the  following  method  finds  the  values  of  the 


we  assume  for  convenience  that  E  is  indexed  =0  that 
its  positive  elements,  lie  along  the  principal  diagonal,, 

K  Select  any  positive  h ;  in  b  for  which  the  subscript 


j  has  not  been  chosen  ore viousl v3  If  none  exist ss  po  to 

instruction  4= 

2„  Locate  the  corresponding  positive  entry  e  ;  1  in  the 

<Xv\dL  *  ve.c  Uv- 

h  marix  and  reduce  the  L  mat  rixf'by  the  Gaussi.  n  reduction 

me' hod  on  the  jili  row  of  li3  Ml  e'ements  of  *  he  jth  row 

become  0  except  the  new  e::,  which  is  1„ 

J  J  1 

30  ttedef  ine  ii  and  >.  to  be  the  m.  t rix  and  vector  result 
'  ing  from  step2„  and  return  to  lc 

A0  The  vtt '  ues  of  'he  r^*  are  read  directly  from  the 
final  h  vec’cr,  D*  is  identified  as  the  sub  form  of  li  whose 
columns  correspond  *o  the  positive  components  of  he  and 
r *  is  the  kth  such  positive  component0  We  obtain  the  cor-- 
responding  lower  bounds  for  w  as  in  lemma  6. 

We  remark  U,a<  in  step  2  above  if  ejj  turns  out  to 
be  nonpositive  the  problem  lacks  finite  optima  and  'he 
t'ru'.Tss  may  be  ?erminatedo 

We  may  illustrate  the  foregoing  method  with  the  same 
noblem  used  to  illustrate  'he  method  of  lemma  5,  We  wr.i‘e 
be  ow  or  ! y  'he  bounding  form,  which  we  have  indexed  \u  cor¬ 
respond  '  o  1  !:e  specifications  of  lemma  7C 
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be  obtain  the  integer  lower  bounds  for  'V^  and  W3  by  round¬ 
ing  the  giver:  values  j">wi  rde  or  =  3,  w-j  -  .5,  In  this 
case  tiie  bounds  are  t-ie  same  as  ft  and  by  'he  method  of  lemma 
-><,  which  we  know  t>a  isfy  he  constraints  of  the  bounding 
fcniJo  :his  will  no1  invariablv  happen,  nor  will  the  amount 
of  corn nu t ut ion  required  by  the  two  me' hods  usually  be  so 
nearly  ’<  he  same*  Generally  sneaking,  the  bounds  implied 
b”  the  method  Gf  lemma  7  v/ill  fall  somewhere  below  those 
implied  by  the  method  of  lemma  5*  On  *he  other  hand,,  f'ere 
mav  be  computational  savings  of  severs 1  orders  of  magnitude 
bv  using  the  bounds  of  lemma  7  to  provide  the  method  of 
lemma  5  with  a  head  start*  fhe  extent  of  the  savings  will 
defend  both  on  fhe  nature  of  the  B  matrix  of  lemma  7  and  its 
associated  h  vector*  In  general,  1  he  method  of  lemma  7 
should  be  bypassed  only  if  h  has  not  more  than  a  single 
positive  component;  hp  and  the  ratio  hj/ejjis  less  than 
or  equal  to  1» 

With  he  fo  lowing  two  definitions  we  complete  r’>e 
inventory  of  the  tools  of  the  bound  escalation  method,  which 
is  presented  immediately  fc» lowing* 

he  will  caH  3  the  dir-.f  ingi-.i  shed  bound"  op  '  om  of  /. 
if  all.  other  bounding  forms  of  A  are  snbforns  of  B0  we  will 
call  1  he  prime  subf-'C m  f  of  A  a  maxima  1  prime  ;nbf ;-rm  if  E 
!n  >  the  same  number  of  col  irons  as  iJ  has  rows,* 


The  Bound  Escalation  Method; 

10  Identity  the  distinguished  bounding  form  B  of  A0 
2U  .'Select  a  maximal  prime  subform  E  of  B  and  apply  the 
procedure  of  lemma  7  to  E,  C The  me+hod  of  lemma  7  may  be 
bypassed  if  h  has  only  one  positive  component  hj,  and  hj/ej j  ~ 
3,  Use  the  lower  bounds  obtained  from  step  2  as  starting 
values  for  the  components  of  x  in  lemma  5,  define  the 
sTarting  value  of  d  ns  d  »  xB„  and  apply  the  method  of 
lemma  5  fo  t Vie  columns  of  B  corresponding  to  E,  If  any 
constraints  of  B  are  left  unsatisfied,  apply  lemma  5  to  all 
of  Bo 

II I  c  Trans  la  t  ion  of  the  Problem  and  Convergence  to  Optima  111- Vo 
We  now  show  how  to  take  advantage  of  the  lower  bounds 
produced  by  the  bound  escalation  procedure  in  the  last  stage 
of  the  algorithm.  Consider  ’he  problem 
(5)  Minimise  wb  +  b^t  •  v/'o 

subject  to  wA  Tp  c  -  w°A,, 

We  relate  (5)  to  (3)  by  means  of 

Lemma  Let  w°  be  any  vector,  and  assume  the  (3)  and  (5) 
have  definite  optima.  Then  if  w  is  a  feasible  (optimal) 
solution  of  (.5),  w*  =  w  «•  vr°  is  a  feasible  (optimal) 
solution  of  (3),  and  conversely, 

Lemma  8  tells  us  that  we  may  let  w°  be  the  vector  of 
lower  bounds  for  w  established  by  the  bound  escalation  pro¬ 
cedure  and  replace  problem  (3)  by  problem  (5)„  Since  w*^  w°, 
this  will  keep  vT  ;  w*  -  w°  nonnegative.  We  know  that 
the  new  c  vector  (equal  to  the  old  c  ~  w°A)  will  be  non- 
positive  for  all  constraints  associated  with  the  old  bounding 
form.  If  anv  components  of  c  are  still  positive  we  may  use 


the  me!  hod  of  l  omnia  [!  to  c re:i ♦ e  a  bounding  form  fha.t  is 

associated  with  a'  least  one  of  these  comnonents,  and  reneat 

rue  process n  Suppose  that  on  some  step  of  exchanging  (3) 

for  (5)  it  turns  out  that  c  -  w°A  becomes  entirely  non- 

positive.  Since  with  each  application  of  lemma  3  we  have 

ataiui  ained  b  0„  and  also  constrained  w  to  be  nonnega*  iva„ 

a  trivial  optimal  solution  to  (5)  is  to  let  0  ~  0a 

1  hus  to  asatue  ilia*  the  algorithm  works  and  that  we  can 

lake  advantage  of  it  if  there  are  two  problems:  we  must  be 

able  i c  force  c  eventually  to  become  nonpositive*  and  we 

must  be  able  t0  salvage  the  w*  which  gives  the  optimal 

A 

solution  fo  the  original  problem  given  w  =  0  in  the  final* 

After  formally  outlining  the  steps  of  the  algorithm  we  will 
address  ourselves  to  these  two  problems; 

General  Form  of  t_he_  Algorithm: 

1,>  Dei  ermine  whether  all  components  of  c  are  nonpositive  <, 
If  so,  we  are  through;  Otherwise, 

2.,  Create  a  bounding  form  by  the  method  of  lemma  4a 
Make  any  legitimate  row  auditions  as  desired* 

3;  Apply  the  bound  escalation  method  to  find  a  lower 
bound  vec*  or  wiy  which  will  satisfy  all  constraints  of  the 
mounding  i’criu 

4,  Redefine  bt,  to  be  w°b  +  l>Qf  redefine  c  ( o  be 
t:  «  w°Ac,  and  return  to  j„ 

kt  j!e  (.V)  nrovi.ks  t  lie  vehicle  for  the  bookkeeping 
of  tt!"  Method,  i all  u  pc:  a  ring  is  handled  by  adding  or 
, ,  ■  b  ;  rae  i  i  ng  ••ows  in  ( 3  ” )  ,  followed  by  subtracting  positive 


integer  linear  combinations  of  the  upper  rows  from  <  lie  bottom 
row., 

1  ho  fo1  lowing  lemma  shows  tha1'  with  the  bookkeeping 
provided  bv  table  (3* )  v  the  problem  of  finding  1  he  optimal 
solu'ion  to  the  original  problem  given  fhat  v)  =  0  in  the 

final  becomes  trivial,. 

Lemma  9,.  '.* he  optimal  solution  w*  to  the  original  problem 
is  the  negative  of  the  vecor  in  the  final  table  in  the 
location  corresponding  to  the  portion  of  the  c  vector  that 
was  originally  the  0  vector  associated  with  the  constraint 
v-i  5*  o,- 

We  complete  the  specification  of  how  convergence  to 
optimality  may  be  guaranteed  with 

Lemma  10*  Assume  that  (3)  has  finite  optima,  and  let  j  be 

a  subscript  for  which  some  c ^  >•  0.  Then  if  A,  bs  and  c  are 

in  ograj  v.  any  rule  for  genera  ting  bounding  forms  which  spec* 

ifa.es  that  the  jth  column  of  A  is  eventually  included  in  a 

bounding  tome  will  assure  "hat  an  optimal  solution  will  be 

found  in  a  finite  number  of  steps*  (The  inclusion  of  c  ^  in 

a  bounding  form  is  of  course  unnecessary  if  c  ,  £ecome$  non» 

J 

poui 1  ive,, ) 

IV  Rules  of  Choice, 

'the  freedom  allowed  for  selection  cixcrr.ativ'c.s 

-it  ■  •  he  algorithm  is  immense,.  We  outline  below  the  various 
provinces  in  which  choice  occurs* 

C-i*  in  the  selection  of  columns  of  A  as  candidates 
for  translation  .in'o  a.  bounding  form,. 


I 


C-2*  In  •'he  se1  ed  5. on  of  elemental  t rans formal  ions 
to  creaie  a  bounding  foriii: 

C-3c  In  the  choice  of  when  to  apply  the  bound 
escalation  met  hod * 

C-4c.  In  the  choice  of  which  maximal  prime  subform 
to  use  in  'he  first  stage  of  the  bound  escalation  method 
when  more  than  one  is  available* 

i'o  obtain  a  problem  solution  expediently  and  efficiently 
there  are  several  considerations  which  suggest  how  the  range 
of  alternatives  may  be  narrowed*  We  examine  the  four  regions 
of  choice  in  more  detail  below,  introducing  such  consider¬ 
ations  as  we  go,- 

C~i  is  perhaps  one  of  the  easier  choices*  Certainly, 
in  the  sf lection  of  columns  of  A  as  candidates  for  irons- 
lotion  .into  the  bounding  form,  we  must  include  one 
associated  with  a  positive  component  of  the  c  vector,  and 
it  is  no1  unreasonable  to  choose  a  set  of  such  columns 
which  are  already  close  to  having  the  bounding  form* 

C-  is  critical*  In  lemma  3  it  was  shown  that  elemental 
transformation?,  of  the  first  type  Crow  subtractions  in  the 
’ abie.-u  matrix)  would  a  1  wavs  suffice  to  keen  the  variables 
normegav  ive*  Hence  elemental  transformations  of  the 
second  • ype  Crow  additions),  which  require  special  cir¬ 
cumstances  before  they  are  permissible,  are  not  strictly 
necessary.,  They  are.,  however,  frequently  desirable* 

Hie  reason  is  as  follows*  he  refer  again  to  the  formulation 
of  Lemma  1  in  which  *  he  transform  R  is  applied  to  problem 
(o)  l o  viel<i  (4)..  Consider  the  effect  of  the  transform  R 


•no  1 

'1  'r-  ;■»  iven  by  7  -  wR  '  w'len  R  is  equal 

fwtec  'r\'J  and  l‘*£  are  inverse  fo  e;  cli 
\  "3  z  :.s  t  :»e  same  as  w  in  every  component 


rs. 

y  When  R  -  1*  "  v.e  have  z 


w 


w 


•  j  " i "if;  c oiipoii.? ri ts  of  z  and  w  are  the  same*, 
i  s  f  i  r -na ;  ton'  o  f  t  h  e  f r  s  t  K  i  nd.  w  e  are 
he  oob1  ron  values  of  i he  variables,, 
of  the  second  kind  we  \ re  generally 
n.sof :•.»  as  raisins?,  *he  solution  values  may 
f  n  is  needed  to  ' oca t e  the  constrained 
dii  ti.cn;;  would  seem  to  be  preferable  'o 
wi  on  •  tey  are  available..,  This  surely  .is 
a  row  addition  i~  is  simultaneously 
r ?rf  cons' rain'  to  he  bounding  form., 
s  .1  o  •  c? xh;  -j  s "  t  h  e  tall  y  o f  c on  side  ra  t  i  <.;  ns 
2.  1. v  ’.ls£-  argument  just  indicated,  it 

two  methods  for  adding  a  column  1  o  the 
haf  reunited  fewer  row  subtractions 
On  tic  o' her  hand,,  tins  is  highly 
•>.  some*  row  rub*: rr c  !  ions  will  create  new 
higher  bounds  ' ban  will  be  created  by 
■jw.lt  •  v  gior'er  sc'  eci  ivity  is  to  t  c\  to 
;  fnr-riice  vtctor  as  large  as  possible  for 
■u:  likely  to  find  their  way  into  a 


.  !,Cii !  a  !  i on  m- thi'vd  is  in  fact  the  pow<  r 
b  .  «-<l  in  he  -re  it:  o*  uer  cons  ieer.v  ions 


should  be  cade  subordina  e  to  getting  as  many  rows  and 


columns  of.  A  into  the  bounding  form  as  possible-, 

C~ 3  is  perhaps  an  easier  choice,  for  which  experience 
may  provide  an  answer,  'the  <",iesri'-r,  is  simply  that  of  wheth¬ 
er  to  apply  ’he  bound  escalation  me*  hod  rs  soon  as  a  bound” 
mg  f>r 'i  is  r.reai  K/J  v  or  *o  vu  i  \  and  try  (o  build  a  larger 
a  rue  tore  for  the  method*  hi  ui'ion  suggests  'hat  for  easy 
problems  Hid  so  n  ion  may  well  be  found  by  a  few  simple 
s  •  ens  bet  ore  a  complex  bounding  fort’  is  created*  On  the 
other  hand,  for  harder  problems  it  couUt  well  be  that  'tying 
•  advance  oward  the  solution  bef  re  a  good  bounding  form 
is  created  would  b?  a  washed  effort  in  'he  sense  tha-  he 
escalation  method  will  take  he  variables  up  1 o  he  same 
pom'  in  roughly  the  varae  amount  of  *  ime  .in  ;  nv  case, 
f,.i.nited  experience  with  problems  a ma  11  enough  to  be  solved 
ty  hand  seems  o  supn-'i*  t  i.ic  notice* 

C-d  ..poear.s  t  tK\s>  print  not  too  critical.  Whenever 
:  vo  or  t-V'ce  iotumr<;:>  of  Che  d is* anguished  bounding  f.-vm  nave 
their  pc;  i*  ive  pompom -its  in  tin  same  row,  one  expedient 
tu  le  fci.  dc ?  ermi  ni.ng  «v.»iicn  one  t  o  use  in  ’he  maximal  s  of  orm 
is  tunpl y  »<>  sclec  •  he  column  assc-cia  ed  wi '  h  the  largest 
*.  1  eoeni  in  'he  c.  vector,,  This  mign<  rt  ?  son;  b  y  be  made 
xubordi;.:  '■  e  to  ?■  r  I’e  which  gives  first  ore1  erence  to  c  -1- 
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.lemma  7  ~o  <  11.  tolnrons  of'  *  he  cli sr  inpuished  bounding  form, 
ra  her  hrn  oiry  '0  .•  'Tine  subform®  Comple  e  reduction 
could  bo  c-  r ?  i «- J  out  using  t.hose  columns  whose  posi'ive 
cnmoouen ‘ s  we "t  unique  to  heir  rows,  until  none  were  lef' 
ha  c-uid  lepimu  el  '  be  used.  i  he  delayed  choice  could 
then  b.  :nat.e  b  rci\:..q'  from  a  se!  of  eligible  columns  hat 
h.-.ve  '  eir  positive  'nnouen  s  in  t  he  same  row  the  column 
so  in  :  ie  *.e  niriology  of  lemra.  7,  =  maxChj/ej- 

where  j  ranges  over  ••  he  columns  indicated,,  (We  reauire, 
of  course*  that  'he  ratio  be  posit i ve„ )  More  refined  rules 
stay  be  invented  if  the  need  for  *  iem  crises® 

Vv  A  Specific  <  igorihm  and  example  Problems.. 

In  order  to  illustrate  4  he  workings  of  the  algorithm 
•/.t  w  i .  i  ;  rbit  tnrily  set  tie  on  a  few  simple  rules  of  chuice® 
i he  specific  method  which  results  mav  be  outlined  as  fo  lows® 
i.  If  the  problem  d-es  not  exhibit  a  bounding  form, 

?.  *>  inmedii;  ;ely  o  instruct  ion  2®  Otherwise p  identify  i  he 
-a.  Xiinal  subioim  H:  i  i'  here  is  a  c  'Oice  to  make  among  more 
ih'.n  uie  column  far  inciui  ion  in  £,  give  firs’  nref  erence 
:  o  t  r  use  not  initially  in  ’he  1  matrix,  and  of  he  remainder, 
sclict  the  one  asscciu  eel  with  he  largest  component  of  he 
vcc •  or®  hake  ;  n v  :  si ? sibte  row  additions  with  respect 

•  i  u 

to  the  co'umns  of  A  associated  with  B®  (If  step  3  has  pre- 
viouslv  he*'!-  curried  out ,  exclude  any  additions  which  would 
r ever  i-  .  subir-ction  performed  in  step  3,}  Apply  the 
b  und  <  id  ion  rou  lore  until  all  constraints  of  the 

o  urmi..1  ;  «  :  r-a  isijed,  and  ad  iu$t  the  C  vector® 


) 


2o  If  c  is  iioiv'fisi ■  ive,  ‘he  problem  is  solved,,  Other- 
vise0  of  i  hose  co:  jmns  j  of  (he  fr  rant  rix  for  which  c  Oj. 
seiect  .he  column  J  which  h3s  the  fewest  number  of  positive 
comoonentSc  If  there  are  ties*  rts1 ric1  J  to  i he  tied  col® 
urns  and  choose  j  so  that  C:  =  taax(c  t)o 

J  J 

Jr.  Consider  oil./  those  rows  of  the  tableau  ma  rix 

for  which  u-^'7'0,  dick  *  he  lexicographically  greatest  row  I 

from  among  them,  arid  for  etch  remaining  row  i  res  ric  ed  es 

; bove  evalu  e  (a)  the  least  (positive  integer)  multiple  of 

row  i  which  when  subtree  ed  from  row  i  will  make  the  result® 

ing  a.,  rionoo.sif  ive„  (b)  the  greatest  multiple  of  row  I 
t  j 

which  can  be  subtracted  from  row  i  and  leave  row  i  lexico» 
graphically  negative  in  the  tableau  ma' rix„  hick  from  (a) 
and  (fc)  the  multiple  that  is  smallest,,  errry  out  the  dncli"tt» 
c.j  sub’ faction  for  fll  .i  as  defined,  and  return  *0  ins' ru8» 

! ion  lo 

/'.ny  ties  not  resc'ved  by  'he  method  may  be  broken  by 
selecting  he  alternative  v.i  h  the  least  index „  Vve  now 
solve  -he  fo1 lowing  example  problems  with  the  method  as  out® 
lined,,  ive  hpve  partitioned  the  tables  to  segregae  he 
ioo  cut  responding  to  the  starting  identity  matrix  in 
order  co  keep  track  of  the  solution  values  of  the  variables,, 
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Co'umii  J  defined  in  s  er>  2  is  indicated  by  the  r rrow„ 
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v-e  .see  s hat  1  r.e  .left  Ik- if  of  ibis  fabt e  is  the  example  pro- 
b  cm  t! ;  <  d  .  o  i  I.  iu.v' rd  e  Hu-  two  methods  underlying  'he 


bound  osca !:•; ioe  procedure-  Identifying  Mu.  ■  wo  by  wo  max¬ 
ima'  :-.'bf  o.rn  in  tiu  .  Ki  rtl  aid  four  ii  columns,  v.u  dread” 

Li  low  t  ha  t  he  lower  bound  j.ou'ion  for  it  is  riven  v:  =  3 

X 

and.  v.„,  :•  5C  . in  _>t  i  ug  he  c  victor  by  sub. riding  ?,  t  imes 

.j, 

the  fir. d  row  a  d  5  imes  ‘he  'hirri  row  we  o  b 1  a  i  n 

5,;  i  .13  -2  -  5  i  ~1  0  -2o 

1  ) 

Mri-ce  re  i  no  -  a  svver  io  ‘he  jroblem  is  vij  =  t„  w-,  -  0C 

,  ZZ  4't  . 

"be  ing  t>i  oblem  was  taken  from'  3  i  ,  where  it 

.  :  u:  cri  ’jv  Gomory  n  cxeur  lify  his  all  in  rger  algori  Ha„ 


■.Us  two  ■iiricre;i  rultj.s  of  choice,  he  solved  ii  in  4 
■i'-'u  ....  nd  3  ■  ivu  s,  rcsoeciivc-.y*  .See  ,fc  -■'•end i x  11  for 
K>r:.  If.  to  comp- i  i  sen  id  f  he  two  m; « iiods0 )  U:  inp  he 

i;is.  <_■•>  of  i;!;j.ico  specified  .  bovc,  wo  ob?  fin  the  s,  me 
vvm  !■  cf  :.b  >  er  ;s  <  lo:uo  r  y  df.;>P  .Mowin'-  fc  r  rcoresoci ;  ~ 

'  ion-  :  d  i  t  Vrt  nf  -v f  o  ig:i  wo  oncia'e  on  them  in  different 
v.i*  -c  >  n  nrosont  -•  ••fob  1  em  for  which  G'-norv’s  m<  hod 
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Ir.s*  e;d  of  using  Che  bound  vscrd  ».tion  method  without  inter- 
ii»"t ioa*  wo  will  work  with  this  t'X"';  tided  (able  -.-Hci  intervene 
v.i.cnf  : -o •  it  !  ir  to  make  legit  its;  e  row  additioriSo  We  then 
obtain  Mir  ro- lowing  sequence  of  tables,  where  at  each  slep 
vie  Sieve  <  or'e - 1 y  r.djus'  ed  i  or  d  vector  according  t  o  the 
previous  fable,,  and  Mien  made  any  legitimate  additions* 
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.  lence  lie  fin;  1  values  for  these  variables  are  (29  30  14), 
which  gives  rhe  id  justed  o  vector  in  table  ft), 

dec.  use  of  the  s*  ace  devoted  'o  their  represent';  t ionB 
i  may  m  cm  at  firs*  glance  ha  >c  r.  ’imu '  a '  i  ons  required 
v.’iiVi  lie  second  stupe  of  he  bound  escalation  method  r  re 
excessive,  dowevrr,  reflec' ion  wi !  show  tha1  a  number  of 
s  eps  of  this  second  s.iage  rn;  v  be  carried  out  in  ’he  amount 
of  time  normally  required  1  or  aj.-io?’  le  pivot,, 

The  preceding  prob’ems  are  n«i  io  suggest  in  any  dog¬ 
matic  fashion  ho  ; he  performance  of  nnr  method  and  Gom'Tv's 
rlgori  hm  arc  til.ely  t  o  coov  re.  A*  Appendix  II  armies, 
ihere  ;  re  reasons  ;o  1 1 1  '  eve  :  hr  '  Gormry ' s  algorithm  mn^ 
pos.i ib  '■  y  be  'ore  efficient  for  simple  problems  which  both 
methods  •*.  n  handle  re!«tiv;-iy  easily,  -iowever ,  some  problems 
which  if:  dif'icul’  :o  handle  wi  h  Gomory's  a’gori’hm  should 
noi  be  too  f'ormid.  ble  with  the  prestnr  code,  ;;  the  List 
examp, e  .how.s,  G»  b,  T'  omps&n  has  found  a  three  variable 
three  in  eon.,  si  y  problem  for  which  Gomorv's  alrorirhm  requires 
over  1  r>0 '■  ’•*ivo‘s.  Our  algorithm  creces  -i  3X3  m:  xima  1  sub- 
form  from  ’he  problem  in  f  nr  s»  eps,  «  no  the  bound  esc;  1 r  t i on 
nr  couf  apri’cd  -q  his  ;.tbform  is  then  able  to  manuf;  cture 
he  so' 1 1)  ion  .in  nc  same  wav  t  ha*  i  handled  1  he  third  ex- 
ampin  pr  >blf'n  above.  In  a  relax'd  vein,  'he  author  lias  made 
;  S'Uiiy  of  Gonorvv;  ..11  integer  •ilr’ori*hm  and  deve’OPOd  a 
si;''o{  ( n  ’a  rw  technique  based  on  concepts  similar  io  those 
un.krl  iiig  '.he  algorithm  of  this  paper  which  significantly 
redicid  ho  number  of  aivo  s  oroirurilv  required  to  ob’.in 
.  ii  oft  :  1.  so  lu  i  ion  ( s  e 


vhis  techninue  m;  v  also 


30 


be  coupled  w’th  t lie  present  algorithm)  as  n  strategy  for 
acceleret ing  convergence. 

He  have,  in  the  foregoing  examples,  used  a  very  limited 
rule  of  choice  in  order  to  keep  the  exposition  simple, 
rn  view  of  lemma  2,,  there  is  a  sense  of  universality  about 
the  elemental  transformations  in  their  applicat't.n  to  this 
problem,  so  that  a  major  avenue  toward  gaining  control  over 
the  integer  programming  problem  may  lie  in  attempts  to 
1(  arn  ttie  best  ways  of  manipulating  there  transf ormations. 
Theorems  about  the  existence  of  bounding  forms,  or  of 
Ijou  tiding  forms  itli  certain  value  for  obtaining  n  problem 
solution,  ant’  the  st-ouence  of  transformations  designed  to 
locate  them,  would  certainly  be  desirable. 

An  alternate  opportunity  for  manipulating  the  problem 
structure  to  obtain  a  good  bounding  form  may  be  to  create 
derivative  inequations  by  adding  together  positive  multiples 
of  i lie  current  constraints.  For  example,,  the  principle  of 
Gaussian  reduction  in  lemma  7  could  be  applied  in  the 
absence  of  a  bounding  form  so  long  as  no  subtractions  or 
divisions  by  negative  numbers  were  permitted,,  Thus  it 
might  be  possible  to  create  a  set  of  constraints  (or  even 
a  single  cons' raint)  considerably  nearer  the  bounding  form 
than  any  in  'he  problem  table,  and  to  use  these  derivative 
constraints  as  well  as  the  original  ones  as  a  guide  for  the 
additions  and  subtractions  of  rows,  and  the  subsequent 
at -era*  ions  of  the  c  vector,  l’he  derivative  constraints 
wou’cl  no  have  to  be  integral,  since  -ihe"  would  necessarily 
corisis*  of  rational  numbers,  and  the  proofs  assuming  integral!  ♦ 
would  remain  valid. 


it  is  conceivable  that  a  technique 


which  combined  these  c  josider:- i  ons  with  a  good  rule  for 
rejecting  iht  e  ■  emt-w  t  ;i  ’  transformations  might  prove  quite 
effective  for  cert;,  in  classes  of  propi  emst. 


af  i' i 


j.e  '.-.if. 


exit  L  r. 


o  r 


!\1Cu 


anr 


(r  5  s),  r]s  ami  l^s  -ire  inverse  to  each  other.-  Since 
each  elemental  transfci~>aiion  consists  entirely  of  integers. 

l 

this  also  imr  lies  thv.i  both  R  and  R  J  aj<  iivepra?,  hence 


'lie  relationship 


w^R  i  gives  that  tit 


r  of  z*  or 


w-v  must  be  .integer  vg  on  .lie  other  one  is*  fbe  direct 
correspondence  of  feasible  solutions  is  observed  immediately 

4 

by  substituting  mr'R " *  In  (4)  under  the  assumption  that  w* 
is  feasible  in  (3),  and  by  substituting  z*R  in  (3)  under 


the  assumption  that  ?*  is  feasible  .m  (4)c  Finally,,  if 


w*  is  optimal  for  (3)*  then  z*  must  be  optimal  for  (-0,, 

•or  else  there  exists  a  z  for  which ‘z(Rb)  +  <  z*('Rb;  +  b  ,, 


A  \  ^ 

and  hence  wb  +  b(.|  ^  w*b  ,  where  w  »  '?,Re  But  by  the 
foregoing  remarks,,  w  is  also  a  feasible  solution  of  (3)„ 
which  is  a  contradiction,,  fl;e  converse  proceeds  similarly... 
r'roof  of  Lemma  2,  We  must  handle  lemma  2  in  several  par  s  . 
Remark  JL  If  R  satisfies  lemma  1  for  ■  he  in'eger  programming 
problem  of  formula i ion  (I)s  for  ail  A°„  bfl  and  c°  satisfying 
■.he  finite  optimality  restrictions,  'hen  R  must  be  integral.-. 


'roof 


and  c 


Prom  formulation  (l )  we  arc  given  that  A  =  (A 
-  (c°  0)„  so  that  we  may  rowri - c  (3)  and  (4) 


i) 


I  v  c  1  y  a  s 

(3a)  Minimize  wb  t  b 

0 

sub.jec*  to  wA°?i  c°  and  w  -y  (),, 
(  ia)  Minimize  ?.  (Rb)  *  b 

0 

subject  in  x(fU°)  >,  cu  and 


.1  l 


iJunpose  that  b.j>  0  tor  each  component  of  b  and  A0  =  R"*„ 

Then  (da)  and  (4a)  become 
(3b)  Minimize  wb  +  b 

subject  i  o  wR“^  ^  c°  and  Y^-  0„ 

(4b)  Minimize  z'JVb)  +  b(J 

subject  to  z Y'  c°  and  zR  2*  0„ 

We  will  show  th3.t  for  the  proper  choice  of  c°  problem  (4b) 
must  have  finite  optima  when  each  element  of  b  is  positive,, 

O  A  .  A  (  # 

Consider  c  =  c  .m  (3a),  where  c  is  integral  and  has 
all  positive  components,  and  let  A°  =  Ic  Then  (3a)  and  (4a) 
reduce  to 

(3c)  Minimize  wb  •-  b,^ 

.  A 

subject  to  w  c  and  w  0, 

(4c)  Minimize  z(Rb)  +  br. 

0 

subject  to  z\X  >  c  and  zR^.0, 

We  see  that  (3c)  .is  trivially  and  uniquely  optimized  for  b.  >0 
by  ie'ting  each  component  of  w  equal  the  corresponding 
component  of  c.  But  then  (4c)  must  have  the  finite  optimum 

A  \ 

given  by  z  =  c  R  Prom  tills  we  may  conclude  that  any 
choice  of  c°  which  leaves  tie  solution  set  of  (4b)  nonempty 
afro  implies  that  (4b)  has  finite  optima*  Evidently  (4b) 
citis !  have  finite  optima  whenever  this  is  implied  oy  the 
con  j'. esc  i  ion  ot  zR '£■ ' '  and  the  objective  function*  3ut 

(4c)  has  the  same  objective  function  and  the  existence  of 

.  .  A 

liiiM  ontima  there  is  implied  by  zR2l  c  and  zR  0o 

hr oce  c  7*  Ov  zR  c  m->i  ies  zR  )?■<>„  and  if  the  latter 
inplus  the  absence  cf  finite  optima  the  former  must  also,. 


an 


But  since  (4c)  does  have  finite  optima,  we  ■■i-tavclude  by 

contradiction  that  (4b)  does  too0 

We  now  show  how  we  may  obtain  a  further  contradiction 

by  (l  S  j  inning  R  is  nonintegral  and  satisfies  lemma  10  If 

R  satisfies  lemma  1,  then  the  optimal  solution  c°  ~  c  R"1  of 

(4c)  must  be  integral  given  that  the  optimal  solution  c  of 

(3c)  is  integrals  We  observe  moreover  that  z  =  c°  «  c'R^ 

is  a  feasible  solution  for  (4b) v  since  substituting  this 

value  in  the  two  constraints  gives  z  -■  c°l?c0t  arid 
A 

zR  =  c  R  -  c  0fc  where  the  inequations  are  more  restricted 
than  we  have  shown  them?  but  iri  ary  case  satisfied,.  Now 
suppose  that  R  is  fractional  in  some  component  in  the  ith 
row  (r.  )0  We  select  c  large  enough  so  that  c  +  (r,  )  >  0, 
and  define  z*  -  c°  +  (0  0oool^o„o0)s  where  as  before 
c°  ~'c  R“*c  Then  is  a  feasible  solution  of  (4b)  since 
clearly  z*  c°„  and  =  c  +  (r^  )  >  0  satisfies  the 

*  o 

second  constraint*  We  note  from  this  last  that  whenever 
c  is  integral?  £*R  is  not„  But  the  feasible  solution  of 
(3b)  corresponding  to  of  (4b)  is  given  by  w*  «  z*R„ 
which  contradicts  the  assumption  that  R  satisfies  lemma  1* 
Therefore?  R  must  be  integral®  (We  note  that  if  c°  ™-  t  A  (Out  • 
then  2*  is  also  an  optimal  solution,  of  (4b)  for  which  the 
corresponding  optimal  solution  of  (3b)  is  nonin tegraL® ) 

Pern ark  2®  R  satisfies  lemma  1  if  and  only  the  determinant 
of  R s,  / R J  -  +1 ,  given  that  R  must  be  integer® 

Proof  If  jR|  -  _+3w  then  R  is  nonsingular  and  has  an  inverse? 
hence  given  that  v  is  integer  lemma  1  must  be  satisfied® 
sin  the  other  hand s  if  K  makes  lemma.  J  true?  must  C/ist, 


XV 


and  by  a  reapplicaf ion  of  the  reasoning  of  Remark  1  we 
know  that  it  must  be  .integral.:,  But  iRj  \R"M  -  jl|  ~  1,. 

Since  the  determinants  of  integer  matrices  must  be  integers^ 
we  have  \R  )  -  |  R"~|  -  +1, 

Remark  3 0  We  need  only  consider  the  case  where  j  ft ,  -  l,, 
since  by  re indexing  two  rows  or  columns  as  the  lemma  permits 
us  to  do,  we  may  change  the  sign  of  the  determinant  if  it 
is  negative., 

Remark  4C  We  will  denote  the  transpose  of  a  matrix  by  the 
prime  (*)  superscript 0  Let  h  be  a  column  vector  of  R  (or 
any  integral  column  vector).,  Fhen  there  exists  a  matrix  X 
which  can  be  expressed  as  a  product  of  elemental  trans¬ 
formations  such  that  h*X’  =  (0  0  ooc  0  k)  where  k  is  positive. 
Proof i  We  need  only  to  show  tha-  the  remark  holds  for 

h“  «  (h  „  h  )0  If  h  and  h  are  both  positive  or  both 
i  2  1  c* 

negative  we  may  proceed  by  aiv&ys  subtracting  the  smaller 

(in  absolute  value)  from  the  larger,;  giving  a  strictly 

monotone  decreasing  sequence  of  absolute  values  until., 

(since  they  are  integer)  one  of  the  components  is  zero. 

If  we  are  left  with  (g  0)„  where  we  do  not  specify  whether 

> 

g  is  positive  or  negative e  we  may  change  it  to  (0  g)  by 
the  sequence  (y  O),  (g  g)u  (0  g)„  using  the  obvious 
additions  and  subtractions.,  If  we  «£nd  up  with  (0  g)  *  (0  *k 
we  may  obtain  (0  k)  by  the  sequence  (0  «k)t,  ( 

( -- k  l-! ) v-  (=k  k)»  (0  If  and  h  arc  of  different 

signs  we  may  first  add  the  larger  in  absolute  value  !o  the 
smaller  and  proceed  as  before. 

Remark  5,  Using  the  method  of  Remark  4V  we  may  create  ar 


X  which  is  the  rmuducr  v  (  elc-monfal  t  ran  si"  onaa  t  i  c-m,  which 
will  i  ran?.;.:  orm  R  into  the  identity  matrix.?  heimna  2  f  o !  3  owe 


i.  .  oiu  ; 

Proof-  Sejt/cr  fn  $\  the  last  column  of  R,  and  reave  e  all 
coEponer. ' s  to  except  *. he  bo i  tom  one  which  we  leave 
'  O8.1  f.  i v c o  flies;  move  iy  the  next  to  the  last  column  of.  ib. 
and  exclude  the  bottom  row f  accomplishing  the  same  resuli. 

•..-•.i  t  h  5.v's  in  all  rows  except  the  next  fo  las'  and  pc  s. sib :  •/ 
the  las'-  Since  rhese  row  additions  and  subtract  ions  do 
not  involve  r he  iast  row,.  none  of  the  0" a  in  the  last 
column  wi  1 '  bo  cl;  an  god  ,,  Repeating  t"is  process,,  eventually 
all  elements  (o  the  right  of  the  main  diagonal  will  be 
zero,  and  all  elements  along  the  main  diagonal  positive 
except  possibly  *  he  one  in  +he  first  row.  be  know  tha  ‘he 
determinant  of  this  final  matrix  R  is  given  bv  }X  It}  •*  jX'flRj 
1»1  3  1.  since  the  de1 eiminant  of  each  elemental  t ran  si onnatio 
i s  .1,.  Finding  the  value  of  [  XRj  by  cofacrors  of  r!  e  last 
column.,  we  see  that  the  bottom  element  on  the  mio  diaronu  ? 
which  is  positive,,  tens  ■  he  1,.  and  its  minor  must  also  be  1  - 


S'ffCeod itig  success!  wi v  i  com  minor  to  minor  v.o  apply  the 
N.ut'.f  /  tguftof-n t  so  sec  tha'  ail  the  diagonal  e ) enen i  $  r"!St 
be  x,  including  finally  the  last...  Huts  v;e  mu  v  tresdi  ly 
n.vKx;  ell  r>  3.!-  ■  ’ling  i/v'i-id i agony  I  on*  ries  in  •  lie  matrix  w vo 
bv  adding  :r  sub  •  rac  i  i.  ag  'he  appropriate  ir.i  eg  or  mul  1 1  pie 


of  el  cctn'r.l  sacs  formations;  cl  such  shat  X  R  *  I,  he*  1  ir.g 
h  »  V  ;  „  t:.r  5 i.  p  -oved 


V3 


Proof  of  Lemma  3,  Since  J’^s  and  f^s  are  inverse  to  each 

other,,  when  R  =  T*s  we  have  z*  =  vtr*  T*5,  and  z*  is  the  same 

as  w*  in  every  component  except  z  =  w  *  +  w  *0  Hence 

S  f  s 

j»2  jl 

z*  is  nonnegative  whenever  w*  is„  When  R  =  z*  =  w*  • 

and  z*  and  w*  are  the  same  except  for  z  *  =  w  *  »  w  * , 

S  o  JT 

Hence  we  must  assure  that  z  *  is  nonnegative  in  some  oilier 

way0  In  (lie  constrain''  z(RA)'&  cB  we  have  the  Jth  column 

of  RA  the  same  as  the  Jth  column  of  A  except  in  the  rrh  row, 

in  which  we  find  the  new  coefficient  (a  .  v  c-  Under 

rj  sj 

condition  (ii)  of  lemma  3  this  coefficient  must  be  nonposa^iv 

as  must  all  other  coefficients  in  the  Jth  column  except 

a  We  may  rewrite  the  constraint  associated  with  column 
5  J 

J  as 


a  z  >  c  +  L(z,  »  k  d  s.) v  where  L  is  a  nonnegative 
s  J  s  "  J  k 

linear  combination  of  the  z,  for  lc  other  than  s„  Since  for 

k 

any  feasible  solution  w*  of  (3)  we  have  z*  =  w  *  for  k  '/  s, 

&  K 

and  since  c.^  0  by  z  *  must  be  nonnegative  in  order 

J  s 

to  satisfy  the  above  constraint,. 

Proof  of  Lemma  40  As  long  as  we  have  two  positive  integral 
elements  in  the  Jth  column  we  may  always  subtract  the  one 
in  the  lexicographically  larger  row  from  +  he  one  in  the 
lexicographically  Waller  row  while  maintaining  both 
lexicographic  ordering,,  and  by  lemma  3,  nonnegativity,. 

Since  ai  least  one  of  the  positive  coefficients  in  the  J*h 
column  is  reduced  by  an  integer  amount  a1  each  step,  as 
long  as  more  than  one  exists,  eventually  all  but  one  must 
become  non pos i  t  .i  vec 


via 


Proof  of  Lemma.  5_n  Because  of  the  form  of  D0  each  individual 

constraint  from  the  constraint  set  w D^d  may  be  written 

in  the  form  d.  ,w.  ■*>  d  ,  *  L(w.  %  k  i)  where  d  is  the 
ij  1  j  k  i  j 

uni  one  positive  coraoonent  in  the  Jtli  column  of  B,  and  L  is 

a  nonnegative  linear  combination  of  the  w  for  k  f  i, 

K. 

Line e  B  contains  exactly  those  rows  of  0  in  which  the 
posit ive  components  appear,,  we  may  rewrite  the  above  as 
b-.x^  dj  +  L^x^,  k  ■/  i)  -  j-  L^Cw^;  /  x^)  where  again 
bj  .  is  the  unique  positive  component  in  the  jth  column  of 
BP  L.  is  a  nonnegative  linear  combination  of  the  remaining 
x^9  and  is  a  nonnegative  linear  combination  of  those 
which  are  not  represented  by  any  x^0  Since  we  require 
any  feasible  solution  to  be  in  integers,,  wc  may  write 

h2'  <(Jj  *  Lo(V  k  '  i)/fci> 

tv c  know  that  each  of  the  must  be  nonnegative,  hence  if 

there  is  any  d.  chat  is  nositive  we  obtain  a  positive 
'  J 

lower  bound  for  the  corresponding  Xj,>  letting  =  0o 

But  as  soon  as  this  lower  bound  for  some  x.  is  known,,  then 

1 

we  may  compute  a  lower  bound  for  the  L  associated  with 
each  of  the  o'her  hence  giving  new  bounds  for  these  x.„ 
as  in  step  3  of  the  method  of  lemma  50  If  we  wish  to 
resist  er  onlv  the  incremental  values  given  io  the  in 

i 

each  step  we  may  redefine  each  d.  to  equal  the  old  cl.  *■ 

J  J  0 

As  long  as  any  of  these  adjusted  d  .  remain  positive  wc 
mav  increment  the  Lower  bounds  of  the  corresponding  x . ,, 

i 

Assuming  the  feasible  solution  set  contains  finite  points, 
the  process  mu  si  eventual  ly  stop  since  i\e  ace  incrementing 


!;c  variables  by  integer  amount s3  At  this  point  the  cons* t  taints 


vii:i 


of  the  original  inequation  wl)  ■?-  d  must  all  be  satisfied,! 

for  otherwise  one  of  the  adjusted  d  .  would  still  be  positive*, 

J 

Proofs  of  Lemma  6  and  Lemma  70  We  combine  the  proofs  of 
these  two  lemmas  since  the  justification  of  the  method  of 
lemma  7  proves  borh0  Vve  assume  that  we  have  adjoined  the 
identity  matrix  above  Es  so  that  when  we  have  completed 
the  Gaussian  reduction  we  will  be  able  to  identify  the 
inverse  matrix  of  the  submatrix  3  of  E  which  has  been 
c’angecl  to  an  identity  matrix0  As  stated  in  lemma  7  we 
assume  that  E  is  indexed  so  that  its  positive  elements  lie 
along  the  main  diagonals  Suppose  h ^  >  0  and  we  carrv  out 
Gaussian  reduction  by  the  method  of  lemma  7  with  the  jth 
column  of  li„  We  will  denote  the  values  of  the  coefficients 
af'er  'he  reduction  by  the  prime  (9)  superscript  Then 
the  formula  for  the  reduction  as  .it  affects  £  and  the 
adjoined  h  vector  .is 


e.  /e . 
ik  jk 


for  lc 


Cl) 


e„,  »  e  "e.  ,/e  . 
ik  jk  ij  jj 


for  k  /  j  ( 2 ') 


where  i  ranges  over  all  the  rows  of  E0  Similarly,,, 

j  (3) 


k 


/h  /e  ,  for  k 
'  k  jk 


n  »  e  °h/e  for  k  f  j 
k  jk  j  jj 


(4) 


Prom  the  fora  of  E,  e 


»  y  0  if  and  only  if  i  ~  kc  We  will 


show  tiiac  .'his  relation  continues  to  hold  after  the  reduction 


with  the  possible  exception  that  e  becomes  nonpositive 

KK 

V 

wren  h  /  U0  Eince  e  is  positive-  obviously  the 
k  "  JJ 


mentioned  relation  is  not  changed  by  (l)c  In  (2)  we  observe 

that  e  is  always  no -posit  ive,  and  e.  .  is  nonpositive 
.jk  iJ 

excepr  wi  en  i  =  j,  l‘hus  e  °  e  /e  .  is  nonnegative  for 

jk  ij  JJ 

i.  r  jo  and  subtracting  it  from  e  leaves  e.  ^  e . ,  fori/ 

i.k  ik 

if  i  =  then  e  .  ~  e  »  e  °  e  /e  .  =■■  0;  Thus  the 

jk  jk  jk  j/  jj 

reduction  insures  that  none  of  the  e.  which  were  originally 

.Ini 

nonpositive  will  later  become  posi’ive*  Since  by  this 

method  we  are  always  dividing  through  some  constraint  by  a 

positive  element  (e..),  followed  by  adding  a  nonnegative 

multiple  (---e ..  )  of  the  resulting  constraint  to  each  cf  the 
jk 

others,  we  arc-  preserving  the  direction  of  the  inequalities 
of  the  constraints  of  each  step*  Moreover,  for  the  same 
two  reasons,  the  inverse  of  the  original  matrix  which  is 
being  implicitly  calculated  in  this  fashion  must  have  all 
normegative  components,  since  the  identity  matrix  begins 
nonnegativec  Finally,  because  the  inequalities  are  preserved, 

I)  8  . 

e  cannot  become  nonpositive  for  b,  >0  or  else  t-Jie  fact 
kk:  k 

that  the  variables  must  be  nonnegative  implies  that  rhe  r  ,K 
feasible  solution  sc  of  the  problem  is  empty,,  Fhus  the 
method  of  lemma  7  is  well-defined* 

The  form  of  the  constraints  implies  as  in  the  proof-.! 
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step 
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!  he  sumo  Cp 

::c  vi.l 

f;d 

'  .  ,  .  .  .  -  . 
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k  c  <  a  i'  i!  v  e  i 
t  ro.'i  (1)  r  1 1 a  t  h 


ie  reduction  step  only  when  i  .  v  c  see 

J 


h  for  k  /  j<.  As  pointed  out  c!  i  .  c 


c 


o  <:  e,  for  k  /  j,  hence  h  /e  >,  h  /c  .  as  c •  ■"lined 
kk  --  kk  k  '  kk  "  k  kk ' 

•Aden  rise  process  is  comp  i  et  ed  wc  a  ;  ).  p.avt  c;-.rr 


..ll! 


he  i 


eu.icti  !>  step  iviU:  every  column  icr  wnich  ' 
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f  i  v  i  ,  so  that,  e  for  these  columns,-  and  i  he  l>  -111*0 


kk 


for  v  1 

K 


simply  piven  ov  y  ">  ( h  /  v 

si  i-  K  / 

i.c  us  denof  e  tie  vector  consist  .inf  of  theca,  positive: 


"k 


bv  d 


;  n-.i 


Ur  el  denote  -he  original  vector  esc  res  'usd  in,' 


:o  !t  s  d  corresponds  o  h\,  By  the  nature  of  the 
baursi.  n  reductin  method,  d  is  the  solution  of  the 
'•••ju.-t  ion  xU  ~  ds  '-.here  3  is  the  submit  trix  of  i:  v.  ,i.ch  is 

j 

v-  us  fe  rnie.-d  into  '  he  identity  matrix  An  the  process  of 
r' -V'”j  up  a  tc  d  ,,  and  x  corresponds  *  0  B  as  v  cor  re  .'■•non-  s 
to  !•»  thus  we  have  d  '  =  d  '*  and  wo  sec  ’hat  is 

;..a  r<-  as  'he  r  vac '  o  r  defined  in  lemma  fcc  nr  have-  a1  read/ 

»  , 

xu.ii  t :;u !  _  ‘  oi  c or: pw.-.e: tv  1=.  Vie 

ris'  f  in.i  1 5  y  shoo  that  no  reduction  wr.iCi  crca'es  an  i  ■  ‘  t. :  •  - 
"  4*  ’  r  i  a  .» 1 1  i*f  a  <’  i  £  t>  reti  r  subtsat  ri  v  of  k  \.i  •  1  imply  a 
nd'.,",  v:i  I  .  1  ■  ''!v  u.-tv  y  , 

lira7,  '.c  ••  :.»«f  rve  tha?  x  =  d  sM  is!  jok  rite  const  r ::  I  r  7 t 

x  :  : ,  T'.'Vibeet  w~  do  not  require  x  '  r.  be  j  a :  *7.  1  a  1  ,  c  1  * 
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/.  A 
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fho/ij  Or.  use  i  :sn  red-  r.  t  ion  rrf  v'cvi'.';  ‘he  direction  of  the 

* 

iiiTOUd  *  i  t  i  es  „  a  nc!  ;:pv  values  for  d  vd  sch  are  obi  ar  n'»{J 
b.'  f  r  a  r :  pg  i 'ip;i ■  i  r  v  matrix  out  of  1  must  f  ive  lepj  tarra 
lo*er  boindr.  !  •  r  ire  aoc'roori s  ■  conpon^n >s  of  if  a>’y 
of  in  I'ali'Or-.  a  a.  rr«  a  !  or  '  iuv--  •  hese  **  ound  b*.'  the  sif  ’  hod 
of  h  'n  nav*.  a  c  .t  *  r.  J  i  r.  t  i on,  since  the  La f  t  or 

v.. !  do  hi  f  act  sa !  i  *  £  v  all  the  c  -nsi  r&ints.  While  seme 
uiV-T  r"  t  beta  <•  r  j.  x  ti  of  -tav  -irmly  ‘b  sane  lower  bounds  as 

hr.;  jdf.r  •  .i 1  joi  b.  t he  ran • '  ••*d  <»f  iemna.  ?,  i  •  cannot  <U> 
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.oo  5.1  ..  i  1  -wn<  s  at  •  art  "ontive  wit hne t  chi?  l.-o ng 
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the  fina;9  t  have  by  lemma  1  that  the  optimal  so  Jut  ion 
<>i  ■  :.c  problem.  pi  ven  by  the  original  tableau  matrix  and 
the  i  i  n  a  I  c  vcC  ui  is  -•  /**<.-  In  this  case,  since  z~“ 
v:*  -■  0  also.,  we  now  seek  a  vector  w°  so  i  hat  the  final, 
c  vector  Cp  and  the  original  c  vector  e  in  a  /  be  related 
bv  c  =  C  «  ,j0Ag  and  the  final  b  eciu&l  t o  riie  or irinul 
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sub  jec 1  *  n  ‘dA  ~a*  c ..  -  w°A 

we  are  assured  that  w  *  *°  is  the  optimal  solution  to  the 
original  problem 

•'ifinize  wb  ♦  b 
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subject  to  wAut  c,„ 
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oince  v.  -  u,  *he  on-isal  solution  for  the  original  profiler; 
i r  r-;r-‘  !y  v.°„  provided  hat  such  a.  w°  exists,  We  riov.  idw.v. 

:hat  this  is  the  case*, 
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'  i  + 1  P 

given  bv  w°  =  w  R  +  W „ R .  *  c.. 0  +  w  R  .  (Hie  identical 
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i'i.i  o t  o:  Ljera-  V\  Hie  ; k r o c ■  f  o!  this  iomr’Ut  .is  given  by 
'J .  iph  (’OW)i\'  o  '.I (-non st fate  the  convergence  of  his  a!  5 
!•!*<  ,-er  i  foci  •  hn  ( sec  j  ?]  ),■  Wu  \*  i  i  i  not  ren/oJ'ir  e  'he 
hr.:1  si?  r  t  he  f  ra as’ at  ions  of  terni oology  fvotor- 
catigsnp  rows  and  columns,  replacing  !  ex  in  op  rnnhic  ;-osi' 
hv  iireai  ivi’V;  sttbs'  ituting  bounding  form  for  pivot  row) 
hu'  straight  ’.forward,  f rora  which  the  applicability  of  the 
proof  as  .it  appears  in  [_ 3j  is  immediate., 
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Gome rv_*_s  A 1 1  Integer  Algora  t Tim ;  Comparison  and  Conj  rasi 
Ait  hough  Gosnorv  develops  his  all  integer  a)  por  i  ■  inn 
in  •<tus  <>  dual  variables  and  pivot  ing  operations,  it  may 
be  explained  very  simply  by  meats  of  the  concepts  present  •.  d 
in  tliis  paper0  Such  an  explanation  will  be  useful  for  a 
clearer  understanding  of  the  relation  of  the  two  algorithms^ 
As  will  be  shown 9  Gomorvts  algorithm  may  be  regarded  essential 
as  a  variant  of  the  present  algorithm  in  which  (i)  freedom 
of  choice  in  applying  the  elemental  transformations  is 
.’-.."••yd  for  a  restricted  routine  to  insure  methodical 
progression  toward  a  solution  with  each  step*  ' afl  ;.?* 
applying  ‘he  restricted  mie„  a  multiple  of  one  of  the 
constraining  inequations  is  added  to  another  to  produce  Tht. 
type  of  constraint  manufactured  by  lemma  4  in  which  exactly 
one  coef 1 icient  is  positive*  (iii)  a* •  row  additions 
possible  in  the  new  constrain'  are  carried  nut,  and  (iv) 

The  lower  bound  implied  for  the  single  variable  with  a 
pc  s it ivc  -oeflicie  t  in  this  now  constraint  is  calculate  1 , 


and  che  c  vector  adjusted  accordingly .. 

.vo  h.  vo  indicated  for  the  firs'  example  p  rob  Ion  in 
ctioi  V  rha4  the  res*  rioted  version  of  cr  method  obtain? 
the  sahi  sequence  of  -  ab‘  es  as  Gomorv’s  method.,  fli « s  is 
vo  out  irol  v  accidental  Pu.  simp'.-  rule  for  tow  sub'  tact  ions 
jeh  ■•■(.  so  lee  led  for  r  example  method  is  the  same  one 
which  Gororv's  me’iixi  use:;  in  *he  first  stage  of  his  >  i  vo  • 
operation.  However-,,  because  the  'wo  inetlnxis  opera'c 

Hat  -'i  i  to  ree  i  '  y  oi  ."’•‘given  "able,  the  actual  row  s  «b- 
i  me  !  l  us  ivai  •  abl  <  *  o  c  adi  :>  av  not  coincide  a  •  •  e r  t  he  i  :  r  s  * 


i : i 

few  s>e:>So  Moreover s  the  rules  for  row  additi  ns  with 
•he  -.vo  methods  are  slightly  different. 

From  c lie  standpoint  of  our  algorithm  these  differences 
are  ■;  i  i na ri  1  y  differences  by  design^  while  for  Gomor/-s 
method  they  are  mainly  differences  of  necessity,  une  reason 
for  this,  as  we  have  observed,  is  the  limitation  Gomory's 
algorithm  imposes  on  the  selection  of  elemental  transformations 
in  the  first  stage  of  the  pivoting  process.  Another  reason 
is  an  equally  strong  limitation  imposed  on  the  second  stage 
of  • be  pivoting  process,  in  which  Gomorv's  method  uses  the 
s'rategy  of  combining  two  of  the  problem  constraints  in*  o 
a  new  one  in  order  to  find  a  lower  bound  for  one  of  -he 
problem  variables*  (id  should  be  emphasized,  in  this  regard,, 
ihat  Gomory’s  algorithm  is  not  actually  divided  into  the 
’’stages"  we  have  identified,  nor  is  if  designed  - o  employ  a 
"strategy"  of  combining  !v.o  problem  constraints  to  find 
lower  bounds  for  certain  variables®  c  have  put  this  con¬ 
s' ructi  u  on  Gomory’s  method  to  explain  it  in  terms  of  he 
ideas  developed  in  ’his  paper,  though  *  he  method 
eve '  v,  ;i  from  a  some  wen*  different  se!  of  notions.:'; 

BoG'i  of  'he  !  imi  tat  ions  we  have  indicate!  rrsup  f  ror- 
concessi  ns  to  the  pivo* inr  rationale  carried  over  from  • he 
Pj:ro!rx  algorithm,  wiiich  requires,  the  c  vector  *n  be  altt”/.: 
a'  ;;ac,i  step.  Outs  ‘he  method  ignores  the  possib.il  iiv  that 
some  synthesis  of  constraints  other  than  the  one  i •  erne  m-s 
may  c-cca  signal  lv  be  more  desirable  for  itig  c  nvt  r  g<  m  e 

'  v;c :  i  a  possibiii  ’  y  is  no  i  only  meaningful  in  our  a  !  gon  '  h  •„ 
but  i:  sugges  re  bv  the  method  of  lemtr..  V„  as  '■  t.  1,  v-  >  r.  .<  • 

It-*  J 


fhc  specific  wav  in  which  Gomorv's  method  acconrv  ishes 
?  he-  second  s’ ape  of  'he  pivot  ing  process  may  be  described 
as  follows..  Hit  first  step  .is  to  select  an  unsatisfied 
constraint  J  and  carry  out  the  row  subtractions  defined  in 
step  3  of  Hie  example  method  in  Section  V„  By  applying 
these  operations  fo  "he  constraint  .  0> which  is  alwavs 
available  implicitly  if  nor  otherwise,  an  inequation  is 
obtained  that  will  have  negative  coefficients  in  every  rov. 
where  the  Jth  constraint  was  originally  positive- -except  0 
of  co-rse,  in  rot.  I0  therefore,  some  ic.-ovej.'.'-lve  multiple 
of  ’his  inequation  when  added  to  the  reduced  jth  constrain' 
will  make  a  new  constraint  with  all  coefficients  nonpositive 
except  in  row  I£  Prom  the  latter  constraint  a  lower  bound 
nay  be  obtained  for  the  net  w  *  and  the  c  vector  adjusted 
accordingly,.  We  present  this  procedure  explicitly  below, 
Goncry 6  s  A 1 1  integer  Algorithm, 

1..  Pick  a  positive  component  c  ^  of  c ,  If  no  more 
exist  the  problem  is  solved, 

2,  Apply  the  row  subtractions  defined  in  step-  3  of 
ike  examp  e  algorithm  of  Section  Vu 

3.  Carry  out  the  transformations  not  only  on  the 

t  .•polar  rablca--  ma>  *-ix,  but  also  on  the  aj  Coined  const  r>- 1 : 
coir  mi  cot  resounding  to  the  inen-atior.  vv  0,  Add  the 
smalles'  nunnegativ;  multiple  of  the  resn ' ting  const  raint 
to  l:e  reduced  Jti.  col  nun  which  will  make  a  tie-,  constraint 


w  ving  exactly  on;  positive  cocfficici 
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run  his  ri  m or*?  clearly  b,s  seen  where  t < •. c*  r  wo 

a  i  •  1. ns  differ-  We  will  now  hazard  a  goc*.*.  or  mo  a  boo  f. 
w-!0  •  nav  hv  c  v  n  cc  •  el  from  this  di  f  f  oj  enc-e..  By  using  a  ere- 
>i*.  t .  j-'im-u  simple  r  ile  roc  reducing  an  unsa  *  isf  ied  const  r.iint 
G<>m  ■>!  ;;  aif-ori!!un  should  be  faster  (provided  we  do  not  use 

at.  equally  simple  ml?)  on  problems  which  both  methods 
find  relatively  cas1  to  solve*  For  harder  problems,  the 
ml.,  five  merits  of  the  ho  methods  will  depend  on  a  number 
o*‘  considerations*  Generally  speaking*  the  efficiency  or 
.i.ier  f  icii  nev  of  our  method  as  compared  with  (Jomory  ‘  s  win 
•tencsul  on  our  abilitv  to  set  up  bounding  forms,  citne-  by 
a  ? log ’  «-wri ndod  strategy  of  arriving  e‘ omental  •  ransfe.  re:- 1  ion 
or,  as  suggest ed  .it:  Section  Vtf  bv  combining  sue'  a.  s,rdo; 
w  •  •  b  r«  •'  echo  i.'iite  f  r  creating  new  const  rain’s-  At  the 


!C bv 

0  1 

i  r,in;t  i 

i  r,  p  ‘ 

he  ma  10 r  part 

of 

cimi 

cc  :  i 

1  b  0 1 

:  ~,.o  •- 

■i'1 

0* 

.r  ne«!; 

sod  W 

"•in  become  <)" 

1 1 

'  S  i  TT) .!. 

Pr 

’0  G-.-t 

a"  no ' 

!.  ,".at 

\\  t 

:  would 

s  cent 

i  ri1  e  )  o  vc1-. 

on 

•  lie  b 

rend 

c  r-  c  .<  1. 

•  rod 

i  n  s  ■  r 

•;i- ; 

:  ef  of. 

>  a  i : ;  i 

up  low.  r  hour 

;dl> 

bv  de 

a  1  i  n  • 

:  i  1  h 

.  1  1*  \  ; 

ri.'.;'  r  . : 

• 

•  irre. 

his  ■■••on  1  d 

be 

an  ••■., 

vast  : 

.I,  .  I' 

•• 

*  ,  •  •if  1 

iX 

1:  )<:•! 

dill;" 

i  o  r:ns  w\!  c.>\ 

Cfr: 

ri ;  e  r 

■  •  ,i  bi: 

'he 

•s  0 

:  x  ■■  7 

rif 

-  c 

cate 

and  then  icier; 

■  if 

v  the 

1)0'.!! 

et :  ;if 

■  rmr 

w  ]  a 

ns  •  it  j 

c  y 

Com: -a cable  di 

Sen ! 

V  ri  r '  '  ; 

re  i  { 

i  r-  i 

r..irc  t 

boss ;  b 

'  e  to 

find  one  con 

s  i  3 

‘  i  nr:  ; 

of  rnn 

:r'e  '  h 

s  \  r.i\ 

1  e  c  •  .■■  ■ 

» ft: 

.  A  i 

Clip  1 

he  same  lines 

y  i 

r  u  i . 

■  i  ■ 
y.*  • 

Will 

'  \  h 

at  the 

f 

vet  dor. 

of  t. 

hoice  avai  ..bl  e 

in  o  • 

r  ~'.'t 

hod 

y  "r-- 

;■>.  Li,-, 

h.l 

lily 

S  V.  1  I 

1  as  an  assi  1 

• 

’ 

i  »  5>  i 

i  ! 

ce  r  •  u 

*  r  • 

r 

'•if: 

• 

*»  !  H  '  w  1  ' 

(c 

n  reach  can  1": 

•:i  i; 

v  •  •  1  0  1  ‘f 

‘  vl  V.* :  1 

J  ‘  V' 

\  ; 

l  \  ■ 

-r- 

.  r  :i 

.ha:  d* 

i  !  o|  -ruble- 

•  <? 

<.  •  ■*  *  • 

'  hr 

1 

tv 


•>l  'll-.’  r. i.'ijmc  nrol'lr-ms  would  Mien  be  delayed  'iiro  t  . 
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